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We provide several applications of a previously introduced 
isomorphism between physical operations acting on two sys- 
tems and entangled states Jl) . We show: (i) how to implement 
(weakly) non-local two qubit unitary operations with a small 
amount of entanglement; (ii) that a known, noisy, non-local 
unitary operation as well as an unknown, noisy, local unitary 
operation can be purified; (iii) how to perform the tomogra- 
phy of arbitrary, unknown, non-local operations; (iv) that a 
set of local unitary operations as well as a set of non-local 
unitary operations can be stored and compressed; (v) how to 
implement probabilistically two-qubit gates for photons. We 
also show how to compress a set of bipartite entangled states 
locally, as well as how to implement certain non-local mea- 
surements using a small amount of entanglement. Finally, we 
generalize some of our results to multiparty systems. 
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I. INTRODUCTION 

In recent years, much of the theoretical effort in quan- 
tum information (QI) theory was focused on establishing 
properties of states and techniques to manipulate them. 
One of the main purposes was — and is — the charac- 
terization and quantification of entanglement properties 
of multiparticle states, as entangled states play an im- 
portant role in several applications of QI. Many schemes 
and applications which involve the manipulation of quan- 
tum states were discovered. Among them, we have tele- 
portation g], purification of noisy entanglement || ||, 
quantum data compression |^| , quantum cloning (?j , and 
quantum tomography Q . 

In practice, entangled states are created by some phys- 
ical action. This suggest that establishing properties of 
operations may play an important role in QI as well. 
First steps in this direction have been recently reported 
[p[dT2[. In particular, in Ref. |l0[ the entanglement ca- 
pability of interaction Hamiltonians between two sys- 
tems has been introduced and analyzed. This quantity 
measures the maximum rate at which entanglement can 
be produced given some particular interaction. On the 
other hand, the entanglement cost for the implementa- 
tion of non-local operations was also considered recently 
]T3| fig] ; in particular, several examples, all dealing with 
an integer number of ebits required for the implementa- 
tion of certain non-local operations, were introduced. 

In Ref. [0 we introduced an isomorphism which relates 



physical operations (completely positive maps (CPM) £) 
and states (positive operators E). This isomorphism 
turns out to be an important ingredient in the under- 
standing of entanglement properties of operations in gen- 
eral. In this paper, we will first review the results ob- 
tained in Ref. [Q: 

• In order to study the separability and entangling 
properties of operations £, it suffices to study the 
separability properties of the associated operators 
E H. In particular, one can use all the results ob- 
tained for the separability of states [0 . This allows 
to answer questions like "Given a CPM £, can it 
be used to create entanglement?". Such questions 
may be relevant in experiments, where one might 
want to know whether a certain machine (set up) 
can be used to create entangled states. 

• One can easily construct physical operations £ 
which can generate bound entangled states (BES). 

• An important problem in the context of distributed 
quantum computation |l8| is the implementation of 
non-local unitary operations. In flj, it was shown 
that an arbitrary 2-qubit unitary operation can 
be implemented using an amount of entanglement 
which is proportional to the entanglement capabil- 
ity of the operation P,[L0| . 

Then, we will discuss several other applications of the 
isomorphism: 

• One can perform two-qubit gates probabilistically 
in the context of single photon experiments via 
creation of entangled states assisted by incomplete 
Bell measurements. 

• Several techniques concerning quantum states — 
e.g. quantum teleportation quantum state 
purification quantum data compression Q 
and quantum cloning — were consider in recent 
years. The isomorphism allows in a simple way to 
obtain similar results for operations. That is, noisy 
unitary operations can be purified, and sets of them 
can be stored and compressed. Furthermore, it is 
possible to clone unitary operations as well as to 
teleport them [ ^9pc| ]. Finally, one can easily see 
how to perform the tomography of general 
non-local operations locally. 

• One can perform certain non-local measurements 
by using a small amount of entanglement. 
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This paper is organized as follows. In Sec. 0, the 
isomorphism between operations and states is reviewed 
and several implications are discussed. This isomorphism 
provides the basic tool for a number of applications which 
are presented in the preceding sections. In Sec III 
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show how to implement non-local two-qubit unitary op- 
erations with unit probability, consuming an amount of 
entanglement which is proportional to the entanglement 
capability of the operation. Sec. IV is concerned with 
purification of noisy operations, while Sec. ^ deals with 



VI 



nt 



tomography of arbitrary non-local operations. In Sec 
it is shown how to implement probabilistically two-qu 
operations in the context of single photon experiments. 
Next, in Sec. VII, storage and compression of no n-local 
unitary operations are discussed, while Sec. VIII is con- 
cerned with the implementation of non-local measure- 
Finally, in Sec 



ments 

to multi-party systems 
Sec. IxL 



IX the isomorphism is extended 
. We summarize our results in 



II. ISOMORPHISM BETWEEN OPERATIONS 
AND STATES 

In p]], an isomorphism which relates physical opera- 
tions [equivalently completely positive maps (CPM) £] 
acting on two systems and (unnormalizcd) states (posi- 
tive operators E) was introduced. This isomorphism is 
an extension of the one introduced by Jamiolkowski [B3[ . 
To be specific, let us consider two spatially separated par- 
ties A and B, each of them possessing several particles^]. 
Let S — {\i)}t—i be an orthonormal basis, and 
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where |<E>) is a maximally entangled state (MES) and P a 
projector on this state. We consider a CPM £ acting on 
the density operator of two d-level systems, one belong- 
ing to A and one to B. Then, there exist an isomorphism 
(linear one to one correspondence) between the CPM £ 
and a positive operator E |jj defined by the following 
relations: 

E Al , 2 , Bl , 2 = £(PAi, 2 ® P Bl , 2 ), (3a) 
SifiAi.Bx) — d 4 tr A23 B 2 _ 3 (E Al: 2 , Bl . 2Pa 3 b 3 Pa2, 3 Pb 2i 3)- ( 3b ) 

These equations have a very simple interpretation: On 
the one hand, Eq. (E3a|) states that given a CPM £, one 



1 We will denote particles belonging to party A by Ai, A2, A3, 
where each of the particles is a ei-level systems. We will also 
use the notation Ai t 2 to refer to particles A lt A2- A similar 
notation is used for particles belonging to party B. 



can always produce the state E associated with £ by ap- 
plying £ to particles A\Bi if they are prepared in the 
state P = Pa 1i2 ® Pb 1i2 - Note that P is a product state 
with respect to parties A and B, while it is a local MES 
in the system belonging to party A and B respectively. 
On the other hand, Eq. (3b) says that given the state E 
(of particles Ai^, -81,2), one can implement the operation 
£ on an arbitrary state p of two d-level systems (particles 
A 3 B 3 ) by measuring the projector P locally in A 2 ^ and 
#2,3 ■ After a successful measurement — where the prob- 
ability of success is given by p = l/d 4 — particles A1B1 
are found in the state £ (p). In summary, a CPM £ can 
be used to prepare a state E, which in turn can be used 
to implement £ with a certain probability of success. 



III. IMPLEMENTATION OF NON-LOCAL 
UNITARY OPERATION WITH UNIT 
PROBABILITY 

In II], it was shown how to implement an arbi- 
trary non-local two-qubit unitary operation with arbi- 
trary high accuracy and unit probability, consuming an 
amount of entanglement which is proportional to the en- 
tangling capability of the operation. Here, we review and 
improve this procedure. To this aim, — as in Ref. 
we consider a family of phase gates 
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where the cr's are Pauli matrices. We show that: 



(i) The operation {/(ajv) can be implemented with prob- 

ability p = 1/2. 

(ii) By applying a finite sequence of operations of the 

form m), each being implemented with probability 
p =1/2 using (i), one can achieve that the opera- 
tion C/(ajv) is applied with probability p = 1. 

(iii) Using gates of the form ([|) with binary angles 
<xn — tt/2 n , one can implement phase gates with 
arbitrary angle a. 

(iv) An arbitrary two-qubit unitary operation can be 
implemented using a sequence of three operations 
of the form U (a) , assisted by local unitary trans- 
formations. 

While (i-iii) are already explained in Q, the imple- 
mentation of (iv) is different to the implementation de- 
scribed in Q. There, an infinite sequence of operations 
of the form U(a) was required in order to implement 
an arbitrary two-qubit operation, while here a finite se- 
quence consisting of three operations suffices. The re- 
quired amount of entanglement is also smaller using the 
new method. 

Since steps (i-iii) will be crucial for the understand- 
ing of some procedures described in later sections, we 
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discuss them again in detail. We start out by show- 
ing (i). First we note that the operator associated with 
the unitary operation C/(ajv) (||) is given by Ea 12 ,b 1:2 = 
\^a N )A lt2 ,B li2 (ip aN \, where 

IVwUi.a.Bi.a = COs(a A r)|$ + ) j4l , 2 |$ + )s 1 , 2 

-ism(a N )\y + ) Al2 \y+) Bl2 , (5) 

and |$+) = (|00)+|11))/V2, |*+) = l®a a |$+) - (|01) + 
|10))/v2 are Bell states. In general, the Bell basis is 
defined as 

|*i 1> i a ) = l®o-ii,< a |* + ). ( 6 ) 

where cr^i = 1, cti.2 = CuC^.i = cr a , 02,2 = o z and 
= l* + ), l*i.2) = l*2,i) = 1*2,2) = 

|$~) are MES. Note that for convenience — to ensure 
a simple notation below as well as in the remaining 
sections — , we use a redundant definition of the Bell ba- 
sis. We consider a situation similar to the one described 
in (^); that is, particles A\ iB\ 2 are prepared in state 
(§), and PA3B3 is the state on which a CPM £ — in our 
case, the unitary operation U(a^) (^) — should be ap- 
plied. Now, a Bell measurement is performed on particles 
^2.3 [-82,3]- Assume that the result associated to l*^,^) 
[| \E r J:L ,j 2 )] is obtained. In this case, the state of system 
A\Bi is proportional to 

£ ® °fujM^(<rti 2 ® of^j) . (7) 

Thus, as a result of the measurement we either imple- 
ment the CPM £, or some unitary operation followed by 
the CPM. We now proceed as follows: In case the result 
of the measurement was ii,«2 [ji 3 J2] ; the local unitary 
operation a il i2 [o~j 1: j 2 ] is applied on A\ \B{\. In case 
that £ is given by the unitary operation t/(ajv) (Q), one 
readily observes that the resulting operation performed 
on pAiB 2 after this procedure will be (i) U(ctN) if i\ = j%; 
(ii) U(aisf)' ! = U(— ajv) if «i 7^ ji- Due to the fact that all 
measurement outcomes are equal probable, we have that 
with probability p = 1/2 the desired operation U{o.n) 
was applied, while with p = 1/2 the operation U(— ccjv) 
was performed, from which (i) follows. 

Before we proceed, we investigate the amount of non — 
local entanglement (between systems A and B) which 
is required to perform the described procedure. The 
amount of entanglement of the state \ip afr ) @ is given 
by its entropy of entanglement 

E(ip aN ) = -x N \og 2 (x N ) - (1 - £jv)l0g 2 (l - X N ), (8) 

where xn — cos 2 (ajv) = cos 2 (n/2 N ). That is, the 
amount of entanglement required to implement the op- 
eration C/(ajv) with probability p = 1/2 is given by m). 
We have that U(n/2) — —icr x ®a x is a local gate and thus 
E(tp ai ) = 0, while E(ip a2 ) = 1, i.e. one ebit of entan- 
glement is required. For N > 2, we have that E(tp aN ) is 



monotonically decreasing with N. The amount of classi- 
cal communication is given by one bit in both directions 
(the value of i\ respectively has to be transmitted). 

Regarding (ii), we have to show how to obtain a prob- 
ability of success p = 1 by making use of the procedure 
described above. Note that with probability p — 1/2, we 
succeed and apply the desired gate, while with p — 1/2 
we fail and apply U{— ajv) instead. Now, if we fail, we 
repeat the procedure but with systems AipBip prepared 
in the state \ip2a N )- With a probability 1/2 we succeed, 
and otherwise we will have applied U{— a at) 3 to the orig- 
inal state instead. We continue in the same vain, that 
is in the fc-th step we use systems Ai^-Bi 2 prepared in 
the state \ip2 k - 1 a N ) so that if we fail altogether we will 
have applied U(— a^) 2 _1 . For k = N we have that 
U{— ajv) 2 _1 = — f/(ajv), and therefore even if we fail 
we will have applied the right gate, so that the proce- 
dure ends. In fact, the iV th step will succeed with p = 1, 
as U(tt/2) is a local gate which can be implemented with 
unit probability and without consuming entanglement. 
That is, a sequence of N operations of the form allows 
to implement the operation J7(ajv) with unit probability, 
which proves (ii) . 

Let us investigate the average amount of entanglement 
which is consumed during this procedure. We have 

E{U(a N )]=J2U) E(^ aN _ k+1 )=a N f N , (9) 
fe=i ^ ' 

where 

1 N 

f N = - V 2 k E(^j ak )</«, = 5.97932. (10) 

TT * — ' 

fe=l 

In (JoJ) , the weight factor pk — {l/2) k ~ 1 gives the proba- 
bility that the fc-th step has to be performed. Thus, we 
obtain E\U{q.n)\ < CKjv/ooi that is, the average amount 
of entanglement is bounded from above by a quantity 
which is proportional to the angle a at and thus — for 
small ctN — proportional to the entangling capability of 
the operation ETJ . The average amount of classical com- 
munication is given by 2 — {1/2) N ~ 2 bits. 

To show (iii), we use the fact that any gate U(a) with 
arbitrary phase a can be approximated with arbitrary 
high accuracy by a sequence of gates of the form U{oin)- 
That is, any angle < a < tt can be written as 

00 

a = ^n fc 2- fc , n k e {0,1}. (11) 

fc=i 

For each k, we have that nu is either "0" — which means 
that the rotation U(ctk) does not have to be performed — 
or "1" — which means that the rotation U(eek) has to be 
performed — . Operations of the form U(a.k) can be im- 
plemented with unit probability using (i-ii). The average 
amount of entanglement consumed to implement U(a) is 
bounded by E < foaOt ebits. 
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Finally, to show (iv), we use the result of Kraus et al. 
There, it was shown that an arbitrary two-qubit 
unitary operation can be written in the following form 



U AB = V ® We~ iH V <g> W, 



where V, W, V, W are local operations and 



H 



E 



5> 

fe=i 



where < p,k < ir/2. We note that 



e -iH =e -iH le -iH 2e -iH 3 



(12) 



(13) 



(14) 



and e~ h are — up to a change of local basis — opera- 
tions of the form (^) for which we already provided a pro- 
tocol [see (i-iii)]. Using this, we obtain that an arbitrary 
two-qubit unitary operation can be performed using a se- 
quence of three operations of the form U (a) , assisted by 
local unitary operations, which proves (iv). The required 
amount of entanglement is bounded by /oo(a«i + + M3) 
ebits. 



IV. PURIFICATION OF NOISY OPERATIONS 

In this Section, we consider purification of a noisy op- 
erations. We will discuss two different scenarios: 

In the first scenario, we consider two spatially sepa- 
rated parties A and B who want to perform a known, 
non-local (entangling) unitary operation U between two 
particles they share. We will assume that A and B are 
only able to perform the operation U in an imperfect 
way. So instead of performing U on their particles, they 
perform some CPM £jj. The problem we pose is the fol- 
lowing: Given several applications of the noisy operation 
£jj and arbitrary local resources, can the parties A and 
B use them to perform the (noiseless) operation U on 
an arbitrary state of two qubits instead? Under which 
circumstances is this possible? In case this is possible, 
we say that the noisy operation is purificable. In Sec. 
1VA we are going to show when and how it is possible 
to achieve this task. 

The second scenario is concerned with the purification 
of an unknown, local noisy unitary operation £jj, where 
we e xplicitly assume a specific form of noise. In Sec. 
EVB, we provide a procedure to implement an unknown 
unitary operation perfectly, given several applications of 
the noisy operation. 

In both cases, it turns out that the isomorphism (JsJ) al- 
lows to use results obtained for purification of states and 
thus for a very simple solution to the problem. Regarding 
the first scenario, the corresponding problem for states is 
the problem of entanglement distillation of mixed states 
[^J. For the second scenario, the corresponding problem 
for states is the purification of a single qubit |2J]. 



A. Purification of a known non-local noisy unitary 
operation 

We consider two parties A and B, who want to perform 
a joint unitary operation U among two particles they 
share. For simplicity, let us assume that U € 5{7(4), i.e. 
the particles are qubits. The parties A and B are only 
capable to perform the operation U in an imperfect way, 
so they perform some CPM £jj instead. For example, a 
noisy iV-qubit operation can be of the form (25) 



£u(j>) = qUprf + ^1, 



(15) 



i.e. with probability q the desired operation is performed, 
while with (1 — q) a completely depolarized state (de- 
scribed by the identity operator 1) is produced. The 
following analysis is not restricted to this specific form of 
noisy operations. 

The operation U is known to both A and B. Fur- 
thermore, they are allowed to use auxiliary systems and 
are able to perform all operations (including two-qubit 
operations) on their individual sites perfectly. In the fol- 
lowing, we are going to show that the noisy, entangling 
operation £u can be purified if and only if the operator pg 
corresponding to £ u (0) is distillable. We also provide 
a practical protocol to achieve this task. 

The purification procedure takes place as follows: 

• £\j is used to create several copies of pg (see ([3a])). 

• With help of entanglement distillation for states, 
out of pf M a number of MES are created. 

• The MES are used to create a set of states of the 
form (@) , either via deterministic state transforma- 
tion (single copy case) or via entanglement dilution 



• Finally, these states are used to implement U with 
unit probability and arbitrary high accuracy as de- 
scribed in Sec. 



III. 



Now we will show that an operation £ jj — where U is 
an entangling operation — is purificable if and only if pg 
is distillable. This can be seen as follows. On the one 
hand, if ps is distillable, one can use the procedure de- 
scribed above to purify the noisy operation £jj. On the 
other hand, if £jj is purificable, this implies that the uni- 
tary operation U can be performed on an arbitrary state 
of two qubits, using a sequence of noisy operations £ u as- 
sisted by local operations and classical communication. 
Since U is an entangling operation, the corresponding 
pure state Ejj is also entangled. That is, the sequence 
of operations £u, assisted by local operations and classi- 
cal communication is capable to create entangled states 
when acting on a certain separable state. Using the iso- 
morphism (H), we can write this sequence of operations 
acting on a separable state in terms of a trace over sev- 
eral operators E% , where local operations in the sequence 
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correspond to separable operators JlJ. That is, the only 
entangled operators which appear in this expression are 
operators ps corresponding to the noisy operation £jj and 
the resulting state is entangled. This implies that from 
several copies of the mixed state ps, an entangled pure 
state can be created. Note that using entanglement dis- 
tillation for pure states || , this implies that one can also 
create a MES. We thus have that ps is distillable, which 
finishes the proof of our statement. 

Since £\j is a general CPM, ps is a mixed state in 
(I 4 (g)(T 4 , where no necessary and sufficient condition for 
distillability is known (see however p7|,p8[ ). It is known 
that non-positive partial transposition of ps is a neces- 
sary condition for distillability, however there are strong 
evidences that this condition is not sufficient |^7| , p8| . Us- 
ing entanglement purification for states, e.g. via the 
methods discussed in p7|-^9|fl , one may be able to obtain 
a MES starting from several copies of ps ■ 

Given the error model (|l^) , one can obtain a necessary 
and sufficient condition for purificability. It turns out 
that unitary operations which are only weakly entangling 
(e.g. operations of the form U(a) with q«1) are much 
more sensitive to noise than unitary operations which 
are strongly entangling, e.g. the controlled-not (CNOT) 
operation [B(J. This means that the tolerable error, spec- 
ified by (1 — q) — such that purification of the noisy oper- 
ation is still possible — in the case of the CNOT is much 
bigger than for U(a) with a -C 1. For the CNOT and the 
error model (|l5|), one obtains q > 1/9 in order that purifi- 
cation of the noisy gate be possible J3l| . For operations 
of the form U(a), one finds q > (16cos(a) sin(a) + 
as a necessary and sufficient condition that gate purifica- 
tion is possible. For a — 7r/2 13 , this value is e.g. given 
by q > 163/164 w 0.994, i.e. less than one percent of 
noise is allowed in this case. 

Note that the process of entanglement distillation in- 
volves two-qubit joint operations as well. The reason 
why we treat these (local) operations differently than 
the (non-local) operation U can be viewed as follows. 
On the one hand, the parties A and B may be spatially 
separated and the interaction between the two parties 
— for example performed through the usage of a (noisy) 
quantum channel — is much more sensitive to noise than 
the local operations performed by only one of the parties. 
On the other hand, each of the parties A and B may be 
considered to possess a single particle only, each particle 
containing several levels. Here, the additional levels are 
used instead of the auxiliary qubits. In this case, the op- 
eration U is concerned with the interaction between two 
different particles, while all local operations (also includ- 
ing multilevel — equivalently multiqubit — operations) 
are operations performed on a single particle, which are 
much easier to implement. For example, using atoms or 
ions with several levels, all local operations can be easily 
performed p3[ . However, controlled interactions between 
two ions/atoms are very difficult to achieve, which leads 
to the fact that two-particle gates are noisy while local 
gates are practically not. Recall that in state purifica- 



tion, it is similarly assumed that local operations can be 
performed perfectly and that it has to be known which 
is the MES which has a large overlap with the mixed 
state the parties share in order that they can distill this 
specific MES. Similarly, knowledge of the perfect unitary 
operation U is required. 

One may also consider that the local operations are 
noisy. In this case, both the process of distillation of 
states and the implementation of the operation U using 
several different states and Bell measurements will give 
rise to some imperfections. The purification of states 
with imperfect means was studied in p^j34| and it was 
found that no MES can be obtained if the local oper- 
ations are noisy and a certain error model is assumed. 
However, one is still able to increase the fidelity, i.e. the 
overlap of the produced state with a MES, where the 
maximal reachable fidelity is determined by the amount 
of noise introduced by the local operations. So, instead of 
producing MES, one produces some mixed state p. This 
state p may then be transformed to a state which is close 
to the ones of the form (||). These states can then be 
used to implement the operation U in an imperfect way, 
since both the states which are used and the operations 
which are performed are noisy. Furthermore, one should 
take into account that a sequence of noisy operations is 
required in order to implement U with unit probability, 
so the errors may accumulate. For almost perfect local 
operations and very noisy non-local operations, one may, 
however, still expect a purification effect. 

B. Purification of an unknown local noisy unitary 
operation 

Here, we consider a party A who wants to perform a 
unitary operation U G SU(2) on a single qubit. The op- 
eration cannot be performed perfectly but is subjected to 
some noise. We will explicitly assume that the imperfect 
operation £\j is of the form ( p^[ ) pq ] with N = 1; that 
is, with probability q the desired unitary operation U is 
performed, while with probability (1 — q) the completely 
depolarized state 1/21 is produced. Here, in contrast 
to the previous discussion, we will assume that the op- 
eration is local and unknown to A (for example, £jj is 
provided to A by a second party via a black box) . Given 
that party A is able to perform all operations perfectly, 
we will show that the unknown noisy unitary operation 
can be purified, i.e. via several applications of £u, the 
noiseless operation U can be implemented on an arbi- 
trary qubit. For simplicity, we assume that the unitary 
operation U is of the form 

U(a)=e- iaa ", (16) 

where a is unknown, however the analysis can be gener- 
alized to arbitrary single-qubit unitary operations. The 
positive operator E corresponding to the imperfect oper- 
ation £{j is given by 
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E = q\* u )(* u \ + 1-^U, 



(17) 



where \$>u) = cos(a)|$ + ) - zsin(a)|* + ). 

We proceed as follows. First, we project E on the 
subspace spanned by {|$ + ), |^ + )} and relabel the basis: 



|o> = |$+), |i) = -<|*- 



(18) 



In case we succeed, which happens with probability (q - 
l)/2, the resulting state will be 



E = \\* u ){*u\ + {1- A)~l 2 , 



(19) 



where \^u) = cos(a)|6) +sin(a)|l) and A = (2q)/(l + q). 
Given N states of the form (19), one can use the proce- 
dure described in pi[ to purify the noisy state, i.e. to 
increase A. For large TV, the average fidelity — that is 
the overlap of the produced states with the state \^>u}— 
scales like Fa 1- 27v3^' whereas the yield — i.e. frac- 
tion of the number of produced states to the number of 
initial states TV — scales like D « A+^i^ @. That is, 
for N — > oo one obtains almost perfect states \if?u) with 
a yield A. Note that the states y&u) are n °t uniformly 
distributed on the whole Bloch's sphere but rather only 
on the equatorial plane. Nevertheless, one can still use 
the same procedure as described in |24j ] , where a uniform 
distribution was assumed. The corresponding values for 
F and D in our case are at least as big as the ones ob- 
tained in [^U , since we have additional knowledge of the 
state, which may be used to further increase F and D. 
Note that in order that purification is possible, we need 
that A > and thus q > 0. So all noisy gates of the form 
( |l5| ) and U given by ( |l6| ) can be purified if q > 0. 

To summarize, we managed to produce an arbitrary 
number of (almost) perfect states \^u) given several ap- 
plications of the noisy operation £jj. Note that \^>u) can 
be transformed deterministically to \^!u) b y un doing the 
basis change (|l8|). From the results of Sec. Ill, we know 
that the state \^u) can be used to implement U with 
probability p = 1/2. What remains is to show that one 
can implement U with probability p = 1. The simplest 
way to see this is the following: If we fail, we try to im- 
plement U again, we make a third attempt and try to 
implement U and so on. Every odd number of steps, say 
2j + 1, we stop the procedure if we have succeeded in 
(j + 1) steps and did not succeed in j steps. In this case, 
we have applied the operation U in total (J + 1) times 
and the operation IP j times, which is equivalent to ap- 
ply the operation U. This is a one-side bounded random 
walk with probability p = 1/2, where one can easily see 
that the total success probability converges to p — 1. Al- 
ternatively, one can also use the operation U to prepare 
states \^u) with coefficient 2 k a, which is possible with 
probability p — 1/2 2 . These states can than be used 
to implement U with p = 1 following the procedure de- 
scribed in M). For a success probability p = 1 — o(e), in 
total o(e _1 )states \^u) with coefficient a are required. 



Alternatively to the procedure described above, one 
may also use a method similar to the one of Sec. |v| 
to implement U given several applications of £u- By a 
sequence of measurements one first determines the state 
E, from which \^u) can be found and used to implement 
the operation U (which is now known to A). 



V. TOMOGRAPHY OF OPERATIONS 

In this Section, we consider the problem of tomogra- 
phy of an arbitrary, unknown non-local CPM. Given sev- 
eral applications of the unknown CPM £ and using the 
isomorphism (||), it is straightforward to completely de- 
termine the non-local CPM by a sequence of local mea- 
surements assisted by classical communication. To this 
aim, we use the operation £ to prepare several copies of 
the associated state E ([5a]). Now, using tomography for 
states ]2lJ , the state E — and thus via (|3b|) also the CPM 
£ — can be determined. 

Next we show that a sequence of local measurements 
assisted by classical communication suffices to completely 
determine a non-local mixed state E (and thus a non- 
local CPM £). Let A and B be two spatially separated 
parties and {A{\ [{Bj}] be an orthonormal |Bq| basis of 



self adjoint operators in A [B\. 
be written as 



Eab — ^ KjAi 



We have that Eab can 



B 



(20) 



where Ay = tr(Ai ® BjEab) is the expectation value of 
the operator A, ® Bj . Now, by measuring the operators 
Ai [Bj] locally in A [B] and using classical communica- 
tion, one can establish the values of all \j and thus the 
state Eab- In case the operation £ acts on two qubits, 
the corresponding state Eab is a state of two four-level 
systems. The set of operators {Ai} can e.g. chosen to be 

Wh 
Eq. 



)0 ~i:i.ii\ — where are defined in Sec. Ill (see 
) — and similarly for {£>i}. 



VI. PROBABILISTIC IMPLEMENTATION 

In this Section, we will show that the possibility to dis- 
tinguish the state |$ + ) from the other three Bell states 
and the capability to produce certain entangled states al- 
lows to implement probabilistically arbitrary two particle 
unitary operations. This has applications in the context 
of single photon experiments, since our method allows to 
implement two photon gates with a certain probability 
of success, which is already sufficient to implement en- 
tanglement distillation. Note that this should be feasible 
even with present day technology 

In the following, we concentrate on two qubit gates. 
Given the results of Sec. [II, one observes that the 



possibility of creating certain entangled states, together 
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with the capability of performing local Bell measure- 
ments allows to implement an arbitrary 2-qubit oper- 
ation That is, the problem to perform two-qubit 
gates is shifted to the problem of 

(i) creating certain entangled states and 

(ii) the capability to perform perfect Bell measurements. 

In the following, we will discuss (i— ii) in the context of 
single photon experiments. 

Regarding (i), in single photon experiments one is al- 
ready able to create certain MES (e.g. via paramet- 
ric down conversion). For example, MES of two qubits 
were created and used in teleportation experiments |57| . 
In addition, the creation of a three qubit Greenberger- 
Horne-Zeilinger (GHZ) state was reported Al- 
though non-linear elements are required in order to pro- 
duce entangled states, it is much easier to use these el- 
ements in such a way that a known state is generated 
rather than using some non-linear elements to perform 
a controlled interaction between arbitrary states. Ap- 
plying the isomorphism (||), one observes that the state 
E corresponding to a general two-qubit unitary opera- 
tion is a pure state of two four level systems (equiva- 
lently of four qubits). For example, the state correspond- 
ing to the CNOT operation [[30| is given by Scnot = 
(|00}A|* + )B + |ll) A |*+) i3 )/v / 2, while the SWAP opera- 
tion (which is given by the mapping \ij) — > \ji)) is spec- 
ified by .©swap = \^ + )a 1 b 2 \^ + )a 2 b 1 - Operations of the 
form (|J) are specified by states of the form (||). Due to 
the fact that the states (JsJ) , as well as Eqnot only have 
two Schmidt coefficients (when considered as a bipartite 
system A — B), it should be possible to create them in 
the laboratory using present day technology. 

What remains is (ii), the problem to perform Bell mea- 
surements. For single photons, using non-linear elements 
only (beam splitters and photo detectors), one is able 
to perform incomplete Bell measurements. In particu- 
lar, one can perfectly distinguish the three sets of states 
{|$+>}, {1*+)}, and {]$"), |*-)} ||. The optimal- 
ly of this process using linear elements was discussed 
in |h|. Due to the fact that Bell measurements can- 
not be performed perfectly with linear elements jlO) (see 
however [^d]]), it follows that two-qubit gates cannot be 
implemented with unit probability using the procedure 
described in Sec. III. However, even incomplete Bell 



measurements (which can already be performed in the 
laboratory) still allow for a probabilistic implementation 
of arbitrary two qubit gates. That is, with a certain 
probability the desired gate is applied, while otherwise a 
different (possibly unknown) operation is performed. In 
the latter case, the input state has to be discarded. 

Let us investigate the consequences of incomplete Bell 
measurements a bit closer. From (||), we know that if 
both parties A and B obtain as a measurement outcome 
the state \& + ), the desired operation was performed. Due 
to the fact that |$ + ) can be perfectly distinguished from 
the other three Bell states using the methods described in 



p9| , p0| , and all measurement outcomes are equal probable 
(in the case of two qubits, p$+ = 1/4), this allows to 
implement the desired unitary operation with probability 
p = 1/16. For unitary operations of the form (^), this 
probability can be further increased to p — 1/4 given the 
fact that also \^ + ) can be perfectly distinguished from 
the other Bell states. That is, if both party A and B 
find as a measurement outcome either \$ + ) or the 
desired unitary operation was performed. In case the 
outcome was additional application of the local 

operation a x is required (see Sec. III). 

Note that probabilistic implementation of two-qubit 
gates is not useful in the context of quantum computa- 
tion, as probabilistic operations may change the complex- 
ity class of the problem and may thus destroy the (expo- 
nential) speed up of the quantum algorithm in question. 
However, probabilistic gates are useful for processes such 
as entanglement distillation Q , which itself is already a 
probabilistic process. For example, this may help in the 
implementation of quantum repeaters |2(| using photons 
only (i.e., for quantum communication over arbitrary dis- 
tances) . Due to the fact that photons are ideal candidates 
for quantum communication (due to their fast propaga- 
tion), it is highly desirable to manipulate them directly 
(e.g. to perform entanglement purification as required 
in the quantum repeater protocol [ ^6| ) rather than map- 
ping their states on the states of another physical system, 
e.g. of an ion or an atom, and vice versa. The method 
discussed in this section may help to achieve this task. 

Recently, an alternative approach was presented by 
Pan et. al p^ ], where entanglement purification with- 
out CNOT operations was discussed. As this approach 
is concerned with a certain distillation procedure only, 
the solution provided in Ref. |^2j to this specific problem 
is more efficient than the one we obtain here. However, 
we provide a more general framework which allows to 
implement arbitrary two-quibit operations probabilisti- 
cally. Another proposal was presented by Knill et. al 
EJ, who showed the implementation of a certain two- 
qubit operation with unit probability, taking full usage of 
all resources (i.e. using an arbitrary number of modes). 

Note that similar techniques may be used to speed up 
— in some sense — slow two-particle interactions. The 
scenario we have in mind is the following: At a certain 
time — e.g. in course of a quantum computation — an en- 
tangling quantum operation should be performed on two 
particles. If the interaction between two particles is weak, 
the required interaction time in order that a entangling 
operation can be performed will be large. Now, instead of 
performing the operation when it is required, we use the 
(slow) interaction at some earlier stage to prepare certain 
entangled states. These states can then — at later time — 
be used to implement the two-particle operation almost 
immediately — once the two particles on which the oper- 
ation should be performed arrive — using the procedure 
described in Sec. III. Although this procedure involves 



local Bell measurements, this will not slow down the pro- 
cess as for the implementation of those measurements, no 
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two-particle interactions are required. For example, we 
can use internal levels of atoms or ions instead of local 
auxiliary qubits (see also Sec. IV). Bell measurements in 



this case involve only single particle interactions between 
the different levels of a particle, which we assumed to be 
much faster than two-particle interactions. 



STORAGE AND COMPRESSION OF 
UNITARY OPERATIONS 



VII. 



In this Section, we will discuss storage [ |43|j20| and com- 
pression of unitary operations. We consider a (possibly 
infinite) set of unitary operations U±, U2, ■ ■ ■ , Upj. Each 
operation is assigned an a priori probability pi. We con- 
sider a long sequence of those operations, where each el- 
ement Ui of this sequence is chosen at random according 
to the probability distribution {pi}. We are interested in 
the average number of qubits which are required to store 
one of the operations Ui and implement the operations at 
later time with unit probability and high accuracy. We 
consider the following variations of this problem: 

(i) The operations Ui are local. That is, a party A locally 

stores a certain number of qubits and uses these 
qubits to implement one of the local operation Ui 
on some unknown state at later time. In this case, 
we are interested in the average number of qubits 
to be stored locally. 

(ii) The operations Ui are non-local. That is, two spa- 

tially separated parties A and B store a set of (pos- 
sibly entangled) states and use these states later 
on to implement the non-local operation Ui. In 
this case, we allow the parties A and B to share 
some initial entangled states. The storage proce- 
dure, however, is restricted to local operations only. 
That is, parties A and B store (and compress) their 
part of the system individually. We are interested 
in the average number of qubits required in A (B) 
to store one of the operations Ui locally. 

(iii) The operations Ui are non-local. In contrast to (ii), 
one of the parties, say A, stores the operations lo- 
cally. Using quantum communication, part of the 
stored system is then transferred to B and finally 
used to implement the non-local operation Ui. In 
this case, we are interested in the average number 
of qubits which have to be stored locally in A as 
well as in the required quantum communication, 
i.e. the average number of qubits which have to be 
transmitted from A to B. 

Note that in all cases, the unitary operation to be per- 
formed is at any stage unknown to A (and B). We will 
show that storage of certain sets of unitary operations is 
possible. Furthermore, the scheme we propose allows to 
compress the amount of required storage qubits (as well 
as the amount of qubits transmitted from A to B in (iii)) 



if one restricts the set of allowed operations to a certain 
subset. It turns out that even for an infinite set of oper- 
ations Ui, the average amount of required storage qubits 
per operation can be much smaller than one. These re- 
sults can be viewed as an extension of the Schumacher 
data compression for states || to unitary operations. In 
fact, we will use the results of Ref. |(| to achieve this 
task. 

Very recently, the problem of storage of a general uni- 
tary operation was considered by Vidal et al. [^0| and 
an optimal solution was provided. In contrary to Ref. 
p(i| ], we propose schemes which are capable to compress 
the required amount of storage qubits and also discuss 
storage of non-local operations. We will propose two dif- 
ferent schemes for storage, one dealing with a possible 
infinite set of unitary operations Ui and one with a fi- 
nite set. We will discuss both schemes in the context of 
(i-iii). 



A. Local storage of local unitary operations 

We start out with (i), the local storage of a set of 
local unitary operations. We consider unitary operations 
acting on two qubits and assume that they are local, i.e., 
both qubits on which the operation should be performed 
are held by the same party, say A. 



1. Storage of an infinite set of unitary operations 

Here, we describe a procedure to store locally a unitary 
operation of the form U(a) (^) with arbitrary, unknown 
a using on average less than four [0.2257] qubits per oper- 
ation if < a < 7r [7r/8]. We assume uniform distribution 
of angles a, i.e. any operation is equally likely. 

We remind the reader that an operation U (a) for arbi- 
trary a can be implemented by a sequence of operations 
of the form U(ctk) @) with binary angles otk = Tr/2 k (see 
Sec. Ill (i-iii)). Using the fact that a can be written 
in binary notation ( pT| ) and assuming that all angles are 
equally likely, it follows that = and — 1 are equal 
likely Vfc. 

We first consider the implementation of U(ak) for a 
certain k = N and — ir/2 N . Following the procedure 
described in Sec. 



Ill (ii), we have that if un = 1, the 



following set of N states is required to implement this 
operation with probability p = 1: 



G N = {\^a N ), |*2a„), • ■ • , l^-OI, 



(21) 



where the corresponding probabilities are given by pi — 
1/2 (_1 for the I th state. If however = 0, no operation 
has to be performed. In this case, one can store the set of 
states Gjv = {l^o), l^o), ■ ■ • , l*o)}, which corresponds to 
the implementation of the identity operation in each step. 
However, each step can be considered independently and 
involves with probability p = 1/2 either the storage of the 
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state \ 1 $ , 2 l a N ) f° r the I th step if njy — 1 or |^ ) if n N = 0. 
Thus, one can us data compression of pure states H for 
each step independently. The corresponding compression 
factor Sj for the I th step is given by the entropy of the 
operator p, which is an equal mixture of the state \^ aj ) 
and l^o)) where j = (N — I). One finds 



Si 



-Xj log 2 (Xj) - (1 - Xj) l0g 2 (l - Xj), 



(22) 



with xj = (1 + cosaj)/2 and Oij = n/2 J . Recall also that 
the Z th step has to be performed only with probability pi . 
That is, the total amount of qubits required to store the 
operation Uiajsi), where it is unknown whether it should 
be performed or not, is given by 



N 



1=1 



N-l 



2>- 



(23) 



We now consider a sequence of operations of the form 
U(a.k) for 1 < fc < oo, i.e. the implementation of U(a) 
with arbitrary a (0 < a < ir). Using (p3]), one finds that 
the total number of qubits needed to store one of those 
operations is on average given by 



(24) 



k=l 



1=0 



k=l 



which can evaluated to be 3.8942. That is, less than 4 
qubits per operation are required on average to store an 
arbitrary, unknown operations of the form U(a). In [ p0[ , 
it was found that on average two qubits suffice to store 
U(a). 

However, if we restrict the possible values of a to 
< a < ir/8 [7r/32], we find that the average amount 
of required storage qubits is given by 0.2257 [0.0206]. 
This can be seen by noting that in this case, the sum in 
Eq. ( pi| ) starts with k = 4 [k = 6] respectively. Thus 
we showed that unitary operations of the form U(a) can 
be stored locally, and that the average amount of qubits 
required for storage can be decreased if one restricts the 
operations to be stored. This result is similar to the one 
obtained by Schumacher [|| for the compression of a set 
of pure states. 



2. Compression of a finite set of unitary operations 

Here, we consider a finite set of unitary operations of 
the form U(a N ) (!), where a N = n/2 N and 1 < TV < M 
and provide an alternative protocol for storage and com- 
pression. This set of operations can be viewed as the ba- 
sic set required to implement arbitrary operations. We 
assume that each of the operation is equally likely. Again, 
we follow the procedure described in Sec. [ni| (i-ii), in 
order to implement a certain operation of the form (|J), 
say U{ajq), with unit probability. The set of states Gn 
( pl| ) is required, where the corresponding probabilities 



are given by pk = 1/2 1 for the fc th state. Note that for 
different AT, a different number of steps are required and 
thus a different number of states has to be stored. As this 
may cause problems, wc fix the number of states to be 
stored for each operation to be M. In case less than M 
steps are required, the state is stored in the remain- 
ing cases, which corresponds to the identity operation. 
Now, the implementation of any operation U(oin) con- 
sists of at most M steps, where in steps (TV + 1), . . . ,M 
the identity operations is performed. The following table 
summarizes the states which are stored for each of the 
operations: 

U( ai ) : Gi = {|* } ), |¥o>, |¥ >, |*o>, ■ ■ ■ , |*o» 
U(a 2 ) : G 2 = {|* f ), |¥ f ), |¥ >, |*o>, • • • , l*o» 
U(a 3 ) : G 3 = {|* f ), |* f ), |* f >, |¥ >, . . . , |* )} (25) 

U(a M ) : G M = {|* ),...,l*f>}- 

Recall that the k th state is always used in the k th step. 
We denote by the fc th column, which consists of the 
k th element of each G; . As the columns Gfc correspond to 
the different steps, we have that column k is only required 
with probability p — l/2 k ~ 1 and all steps — and thus all 
columns G& — can be treated independently. That is, 
we store each of the columns Gfc independently. Due to 
the fact that all states within each column Gfc are equal 
likely and non-orthogonal, one can use data compression 
Q. The compression factor Sk for column Gfc is given 
by the entropy of the density operator pk , where pk is an 
equal mixture of all the states of column k, and column 
k is only required with probability p = l/2 k ~ 1 . Thus the 
total number of qubits required to store one of the opera- 
tions U(a N ), < N < M is given by Y^ =1 S k /2 k ~ 1 . For 
example, for M=100 [1000] we obtain an average amount 
of 0.245 [0.0361] qubits which have to be stored on aver- 
age to implement one of the 100 [1000] operations picked 
at random. 



B. Storage of non-local unitary operations in A and 

B 

Here, we consider (ii), the storage of a set of non-local 
unitary operations. We will discuss variations of both 
protocols described in the previous section, taking into 
account that we now have two spatially separated parties 
and the operations are non-local. That is, the states 
to be stored are entangled states and we consider local 
storage of the subsystem belonging to A (B). This means 
that both, the coding and decoding procedure has to be 
local, but may be assisted by classical communication. 

We first consider the storage of an i nfinite se t of unitary 
operations of the form U(a) (see Sec. VII A 1 ) . W e follow 

however 



VII A 1 



the same protocol as described in Sec. 
we now use a different ki nd of d ata compression. The 
protocol described in Sec. VII A 1 involves storage of two 
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equal likely states, |\I>o) or for some (Xn = tt/2 n . 

Note that the state l^aw) 1S an entangled state, so in 
contrast to the previous section, we cannot use standard 
data compression for pure states, as we are restricted to 
local operations only. However — as shown in Appendix 
A — it is also possible to achieve local data compression 
for a set of entangled states. That is, each of the parties 
manipulates only its own subsystem and can thereby re- 
duce the average number of qubits required to store its 
part of the entangled state, without affecting the entan- 
glement with the other system. Note that this problem is 
equivalent to the data compression of mixed states with 
commuting density operators, where the entanglement 
with some other system should be preserved. It turns 
out that the compression factor for A (B) is given by the 
entropy of an operator p, which is an equal mixture of 
the reduced density operators p\ (p^) corresponding to 
the states l^o), Y&u N )- Note that this corresponds to the 
upper bound on the number of qubits to be stored in case 
the entanglement with another system is not required to 
be preserved H|. While it is known that this is not the 
optimal compression rate in case the entanglement with 
some other system is not required to be preserved, it is 
not clear whether the compression rate is already opti- 
mal under this stronger restriction. In our specific case, 
we obtain 



S 



N 



-x N \og 2 (x N ) - (1 - x N ) log 2 (l - x N ), (26) 



with xn = (1 + cos 2 ajv)/2. Using now this local com- 
pression protocol inste ad of Sc humacher's for pure states 
in the protocol of Sec. VII A 1 , one finds that the the av- 
erage number of qubits which have to be stored locally in 
A (B) is given by 4.7758 if < a < ir. If we restrict the 
possible values ofato0<a<7r/8 [7r/32], we find that 
the average amount of required storage qubits is reduced 
to 0.3976 [0.0379]. 

Regarding the storage of a finite set of unitary oper- 



ations (see Sec. VII A 2), we follow the same protocol 



as described in Sec. VII A 2 and use again a different 



kind of data compression due to the fact that we are re- 
stricted to local operations. This time, data compression 
for a finite set {| V F;)} of M entangled states is required. 
The entangled states are all equal likely and are of the 
form I^q.^). It turns out (see Appendix A) that one can 
achieve a compression rate which is given by the entropy 
of a density operator p, which is defined as an equal mix- 
ture of the reduced density operators p\ corresponding 
to the state One finds that the total number of 

qubits required on average to store one of the operations 
U(a N ), < N < M locally in A (B) is given by 0.333 
[0.050] qubits for M=100 [1000]. 



C. Storage of non-local unitary operations in A 

Finally, we consider (iii) the local storage of a non- 
local unitary operation in A. That is, we consider a local 



memory (in A only), but we want to implement the oper- 
ation non-locally. It turns out that this problem is a triv- 
ial combination of the previous two problems. We ha ve 
that one can use the methods described in Sec. VII A to 



store the operations locally in A, and one obtains the the 
same compression rates. The average amount of quan- 
tum communication from A to B — which is required to 
implement the operation non-locally — can be found us- 
ing the method described in Sec. VII B| . That is, one part 
of the entangled system is compressed and send through 
a quantum channel to B. The compress ion ra te can be 
calculated in a similar way as in Sec. VII B| , however 



one has to take into account that the state |vE , 7r /2) is a 
separable state and thus no quantum communication is 
required to transmit one part of this state. For exam- 
ple, one finds in case of an infinite set of operations of 
the form U(a) with < a < ir [tt/8] that the required 
amount of quantum communication from A to B is given 
by 2.7758 [0.3976] qubits. 

This last method clearly distinguishes between the re- 
quired amount of local storage qubits and the non-local 
content of the operation, i.e. the average amount of quan- 
tum communication. Note that storing the operations lo- 
cally (see Sec. VII A) requires a smaller amount of stor- 



age qubits than st oring a non-local operation directly in 
A and B (see Sec. |yiI B|) . 



VIII. NON-LOCAL MEASUREMENTS 

In this Section, we consider the implementation of non- 
local measurements. We consider two spatially separated 
parties, A and B, each possessing a <i-level system. The 
two parties want to perform a complete, joint measure- 
ment on their system, specified by a set of rank rik Projec- 
tors {P k } such that Y,k P k B = 1a ® Ib< The questions 
we pose are the following: How can the parties implement 
this non-local measurement? What are the entanglement 
properties of those measurements; that is, (i) what is the 
amount of entanglement required to implement a certain 
measurement? (ii) what is the average amount of entan- 
glement which can be produced given a single application 
of the non-local measurement? 

We provide several procedures to implement arbitrary 
non-local von Neumann measurements and discuss their 
entanglement properties. We show that the required 
amount of entanglement depends on the measurement 
to be implemented. We introduce examples of non- 
local measurements which can be implemented using less 
than one cbit of entanglement. One can easily generalize 
some of our results to implement also arbitrary measure- 
ment, described by a positive operator valued measure 
(POVM), i.e. a set of positive operators 0^ B such that 

First, we note that the amount of entanglement re- 
quired to implement the non-local measurement depends 
on (i) whether one is only interested in the measurement 
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outcome or (ii) the system should in addition be in the 
corresponding state after the measurement. For exam- 
ple, one can perform a complete Bell measurement (i.e. 
a measurement in the basis (||)) on a state of two qubits 
using one ebit of entanglement regarding (i), while two 
ebits are required in case of (ii). 

Proposal 1: A trivial procedure to perform an arbi- 
trary bipartite measurement is the following: The state 
of system B is teleported to A, consuming log 2 (d) ebits. 
Then, the measurement is performed locally in A, which 
already suffices in case of (i). Regarding (ii), one also 
has to teleport the particle back to B, consuming again 
log 2 (d) ebits. Note that in case of a complete Bell mea- 
surement, i.e. a measurement in the basis (||), where each 
basis state is a MES, this procedure is in fact optimal. On 
the one hand, one consumes two ebits to implement the 
measurement. On the other hand, one can also obtain an 
average amount of two ebits given a single application of 
a non-local Bell measurement. One just has to consider 
the operator 2% (||) associated to each possible outcome 
of the Bell measurement. One observes that the non- 
local entanglement of all E^ is given by two ebits, and 
each measurement outcome is equal likely. This leads 
to an average amount of entanglement of two ebits. We 
have that if the amount of entanglement required to im- 
plement an operation, E- im equals the amount of entan- 
glement which can be obtained given a single application 
of the operation, E CT , the first process is optimal and 
E{ m is the minimal amount of entanglement required to 
implement the operation. This is due to the fact that 
E CI < -Eim, otherwise one could create entanglement for 
free. However, if one wants to measure the joint sys- 
tem A,B in a basis which is not maximally entangled, 
one might expect that the required amount of entangle- 
ment is smaller than 21og 2 (c2) ebits. With the following 
method, we show that this is indeed the case. 

Proposal 2: We consider the situation were all Pk 
are rank one, i.e. = 1 and thus Pk = \4>k) AB(4>k\- We 
define a non-local unitary operation U by 



(i), the procedure ends at this point. Regarding (ii), A 
and B in addition implement the operation U' to ensure 
that the system is also in the required state after the 
measurement. Alternatively, they could also prepare the 
measured system in state \4>k)j as due to the implementa- 
tion of the measurement, any possible entanglement with 
some auxiliary system is destroyed anyway. We note that 
the choice of the local basis in A, {|o,-)} and B, {\bi)} is 
not fixed and may also change the entanglement proper- 
ties of the operation U. This can be seen by considering 
the following trivial example: We have d — 2 and \4>oo) = 
|00>, |0oi> = |Ol),|0io) = |10) and |^ n ) = |11). By 
choosing |oi)|&i) = |0) and |a 2 )|6 2 ) = |1), we have that 
U = Iab, i.e. no entanglement is required to perform 
the measurement. If we however choose the mapping 
|0oo> |OO),|0oi> -> |1O),|0io> |Ol),|0n) -> |11), 



we find that the operation U = U s , 



which requires 



two ebits to implement pl|. In this case, the choice of 
the proper local basis is trivial, however we do not know 
the optimal choice for a general measurement. Note also 
that this procedure fails to implement non-local measure- 
ments where the rank of some projector Pk is larger than 
one. For example, if Pi = 100)^5(001 and P 2 = 1 AB -Pi, 
this procedure fails to project in the subspace spanned 
by P 2 , as it already gives a fine-graining within this sub- 
space, which is a different problem. The next method 
will overcome this limitation. 

Proposal 3: Here, we consider a complete set of M 
non-local projectors Pk which might have arbitrary rank 
nk- Clearly, ^2kLi n k = d 2 . Alice uses an M-level auxil- 
iary system initially prepared in state |1), which is used 
to label all possible measurement outcomes. We define 
a unitary operation U acting on the auxiliary system A 
and the joint system AB as follows: 



M 



(i| + |i)i01)®P 



MS 



\j)A(j\)®Pf B ] 



(28) 



AB \ 



(27) 



fe=i 



where \k) A B = \ak)A\bk)B, and {\di)} [{\bi}}\ is some lo- 
cal basis in A [B] respectively, with 1 < k < d. The 
procedure takes place as follows: First, the parties apply 
the non-local unitary operation U, using e.g. the pro- 



cedure described in Sec. Ill for d = 2, consuming an 



amount of entanglement which is specified by the oper- 
ation U. In case U is only weakly entangling, e.g. if 
(k\<f>k) ~ 1 (i.e. \4>k) are only weakly entangled states), 
the required entanglement is small (see Sec. |||) 
Then, parties A and B both perform a local measure- 
ment specified by projectors on the states [{|&t)}] 
respectively and communicate the outcome of the mea- 
surement classically. In case they obtained the outcome 
ak , bk , they know that the outcome of the measurement 
is k, i.e. they measured the projector Pk- Concerning 



After application of U , the auxiliary system A is mea- 
sured in the basis If the outcome k is found, one 
readily observes that this corresponds to measuring the 
projector Pk on the system AB. Note that no further op- 
erations are required, as the system AB is already in the 
appropriate state (ii). The amount of entanglement re- 
quired to implement the non-local measurement is again 
specified by the operation U . 

For example, if d = 2 and P x = (|00)(00| + |H)(H|), 
P 2 = (|01)(01| + 1 10) <10 1 ) , it turns out that one can cre- 
ate one ebit given a single measurement of this kind. 
To see this, we prepare the system AB in the prod- 
uct state p = 1/2 (| $+)($+ 1 + |#+)(*+|) and perform 
the measurement. In case we obtain outcome "1" ["2"], 
the state after the measurement is |$ + ) [|^ + )] respec- 
tively. In both cases, we created one ebit. However, it 
is not clear whether one ebit of entanglement also suf- 
fices to implement the corresponding unitary operations 
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U = 1 A ® Pf B + af® P 2 . Although the state E v 
associated to U via (g) has an amount of entanglement 
of one ebit and U = W , it is not clear whether a single 
copy of the state Exj suffices to implement U. 

It would be interesting to establish the minimal 
amount of entanglement required to implement a gen- 
eral, non-local measurement. 



IX. MULTIPARTY OPERATIONS 

In this Section, we generalize some of the previous re- 
sults to multiparty systems. We consider several spatially 
separated systems A, B, . . . , Z, each possessing several 
d-level systems. We first generalize the isomorphism (||) 
between CPM £ and positive operators E to multiparty 
systems. Here, £ acts on several d-level systems, one 
located in each site A, B, . . . , Z and E is a positive oper- 
ator on the Hilbert space Ha 1i2 ® ■ ■ ■ ® Wz 12 . We have 
that Ti-Ai =d d and similar for the remaining parties. For 
a TV-party system, it is easy to show that 

E A ia ...z l - 2 = £(Pa 1j2 <8> • . . ® Pz 1:2 ), (29a) 
£{PA x ...Zi) = d 2N tr A23 ... Z23 

(Ea 1 , 2 ...z 1i2 Pa 3 ...z 3 Pa 2 ,3 ■■■Pz 2>3 )- (29b) 

The interp retat ion is similar to the one of Eq. (^). On 
one hand, ( [29a|) states that E can be created from a TV- 
party product state, wher e eac h party prepares locally a 
MES. On the other hand, (29b) tells us that given E (par- 
ticles A\^,P\,2 ■ ■ ■ -^1,2), one can implement the multi — 
particle operation £ on an arbitrary state p of N d— level 
systems (particles A3B3 . . . Z3) by measuring locally the 
projector P (||) on particles ^2,3,-82,3, . . . , Z 2 ^ in each 
of the locations. Note that the probability of success is 
given by p — l/d 2N . 

As in the bipartite case, one may ask for a certain map 
£ whether it is capable to create entanglement. Since for 
multiparty systems, there exist many different kinds of 
entanglement (see e.g. p6| -|49|] ) , one may also ask which 
kind of entanglement can be produced. Again, all these 
questions can be answered by establishing the entan- 
glement properties of the operator E associated to the 
CPM £ via (29a). In particular, if E is bound entan- 
gled f48|| , then £ can only create BES. In a similar way, 
given some BES one can easily construct the correspond- 
ing map which is capable of generating BES of the same 
kind. 

One may also consider the implementation of arbitrary 
iV-qubit operations with unit probability. On one hand, 
any TV-qubit operation can be written as a sequence of 
bipartite CNOT operations and single qubit unitary op- 
erations, for which we already established a protocol. On 
the other hand, we may consider iV-qubit unitary oper- 
ations of a specific form and show directly how to im- 
plement them with unit probability given certain states. 
We consider a unitary operation of the form 



Un(um) = e 



(30) 



where <Xm — 7r/2 M . It tu rns out that a natural extension 
of the protocol of Sec. [ill] (i-iii) allows to implement 
operations of the form fl30|) with probability p = 1 . The 
operator associated with the unitary operation Un{oim) 
is given by E Al _ 2 ,...,z li2 = IVwXVwL where 

IVw) =cos(a M )|$ + ) j4l , 2 |$ + )i3 1 , 2 ...|$+)z 1 , 2 

-isin(a^)|* + )A lja |* + )Bx, s ■ ■ ■ \^ + ) Zl . 2 - (31) 

Regarding (i), we just note that Bell measurements and 
the corresponding local unitary operations are performed 
at all location A, B, . . , , Z. For all possible measurement 
outcomes, it is easy to observe that the operation per- 
formed on some state p Al ,...,z 1 will either be (i) [/(ckm) 
or (ii) U (-qm), each possibility appearing with probabil- 
ity p — 1/2. Steps (ii) and (iii) can be adopted without 
changes, which finally allows to implement an operation 
of the form (|30| ) with arbitrary angle a and unit prob- 
ability. Note that operations ( |30| ) are capable to create 
GHZ-like entanglement and are thus truly iV-qubit en- 
tangling operations. 



X. SUMMARY 

To summarize, we have provided several applications 
of an previously introduced isomorphism between opera- 
tions an states. First, we discussed how to use this iso- 
morphism to establish separability and entangling prop- 
erties of operations £ and to construct physical opera- 
tions which are capable to create bound entangled states. 
In addition, we showed how to implement an arbitrary 
non-local two-qubit operation consuming an amount of 
entanglement which is proportional to the entangling ca- 
pability of the operation. 

Then, we have shown how to implement several tech- 
niques developed for states — such as purification or data 
compression — also for operations. In particular, we have 
shown that a known, noisy, non-local unitary operation 
as well as an unknown, noisy, local unitary operation can 
be purified. In a similar way, we use these results to 
establish tomography of arbitrary operations. Then, we 
showed that unitary operations can be stored locally and 
non-locally and that the amount of required qubits for 
storage can be decreased, which can be viewed as gener- 
alization of data compression to unitary operations. In 
this context, we also provided a protocol which allows for 
local data compression of a set of entangled states. Note 
that it is straightforward to obtain a number of other re- 
sults which were developed for states also for operations. 
For example, it is easy to show that also unitary oper- 
ations can be cloned (via cloning of the corresponding 
state E) or teleported (via teleportation of the states re- 
quired to store the operation) [19pC|]. In case of cloning, 
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one has to take into account that the cloned states al- 
low for a probabilistic, imperfect implementation of the 
required operation only. 

We also provided a method to implement arbitrary 
two-photon gates probabilistically with present day tech- 
nology, which opens the way for practical quantum com- 
munication over arbitrary distances. Finally, we dis- 
cussed the implementation of non-local measurements 
and generalized some of our results to multi-party sys- 
tems. 
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APPENDIX A: LOCAL DATA COMPRESSION 
FOR A SET OF ENTANGLED STATES 

In this appendix, we consider the problem of local data 
compression of a set of pure, entangled states, where all 
reduced density operators commute. Note that this prob- 
lem is equivalent to the problem of data compression of 
a set of commuting mixed states under the restriction 
that entanglement with some other systems should be 
preserved. Let G = be a set of L pure states, 

where 

|*i) = c ai |00) AB +s Qj |ll) AB . (32) 

and c Qi = cos(ai), s ai = sin(ai). Each state is assigned 
a prior probability pi. Two spatially separated parties 
A and B are fed an unending sequence of states 
where each successive state is chosen randomly and in- 
dependently from the set G according to the probability 
distribution {pi}. A sequence of length N is of the form 

|^ li3 ... iw ) = |* il )|^ 2 )...|*«), (33) 

and appears with probability p ili2 ...i N = PhPi 2 ■■■Pi N - 
The parties A and B want to store the sequences locally, 
i.e they are allowed to perform local operations and clas- 
sical communication. We are interested in the average 
amount of qubits per signal state which are required in 
A (B) to store the signals faithfully. We will use as a cri- 
terium the so-called GLOBAL-FID criterium Q; that 
is we require that the average global fidelity of all pos- 
sible sequences is 1 — e. Note that we consider the so 
called "blind case" [p4[, that is neither A nor B know 
the specific sequence (|33|). 



Let p A = trs(|^E'i)(\E , i|) be the reduced density opera- 
tor of system A of the state \^i) and 

~p A = Y d PiPt (34) 
i=i 

be the weighted average of the reduced density oper- 
ators of our signal source. We denote by S(p A ) — 
tr(p A log 2 p A ) be the von Neumann entropy of p A . 

Given a sequence of length N, N sufficiently large, we 
provide a protocol with the following properties: 

(i) The required amount of storage qubits in A (B) is 

given by NS{p A ) + 5. 

(ii) The average global fidelity (averaged over all possible 

sequences) F is given by 1 — e. 

We have that S is some function which is of the form 5 = 
fiN 13 for some // > 1, 1/2 < j3 < 1 and e — > as N — > oo. 
That is, on average S(p A ) qubits per signal state have to 
be stored locally in A (B). Note that we do not claim 
that this the optimal compression rate achievable. 

For pedagogical reasons, we will proof our statement 
in the simplest case, where the set G consists of two pure 
states only. We will even assume that |^i) = |00) and 
|*2> = c CK |00)+s ce |ll) and pi = p 2 = 1/2. Note that the 
proof can be easily generalized to an arbitrary number of 
signal states and an arbitrary probability distribution. 

We have that 

/^i±£l| )(0| + i^|l)(l|. (35) 

We define local projectors Pa (Pb) acting on N qubits 
as follows: 

fc+ 

Pa = Pb = J2 ( 36 ) 
k=k- 

where fc± = (1 + c 2 a )/2 ± pN? , p > 0, 1/2 < (3 < 1 and 

Pk= |0)(0|® fc ® 11X11®^. (37) 

perm 

The sum in (j3?]) runs over all possible b^^ = N\/[k\(N — 
k)\] permutations (without repetitions) of k zeros and 
N — k ones. Thus P^ is a projector in the subspace 
spanned by all states which contain exactly k zeros and 
(N — k) ones. The dimension of Pk is given by 6jv,fe- 

The projector Pa (Pb) is measured locally in A (B). 
If the measurement is successful, \og 2 (d) — where d = 
dim(P J 4)) — qubits are used to store the resulting state 
in A. This can be accomplished by relabeling the states 
which span Pa to {|Z)}f=i an d store those states lo- 
cally, which clearly requires log 2 (d) qubits. The decod- 
ing procedure consists of undoing the relabeling. In case 
the measurement is not successful, some state \0e) is 
stored instead. We show that: (i) log 2 (d) = NS(p) + S 
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and (ii) F 



T,i 1 i 2 ...i N Piii2-iN F i 1 i2...i N > 1-c, where 

Hi2-m = \(^iih-m\ P A ® -PB|*iii 2 ...iiv)| 2 and tne sum 
runs over all possible sequences. That is, the storage pro- 
cedure requires the announced amount of qubits and the 
average fidelity is sufficiently large. 

Regarding (i), it is easy to see that d = dim(P,i) < 
(k~ + k + + l)6jv,/t , where ko £ [k~, k + ] is the value that 
maximizes 6jv,/c in this interval. Substituting the values of 
k~ , k + in this bound, we find that \og 2 {d) = NS(p A ) +5 
as required. 

We now concentrate on (ii) , the average global fidelity 
F. Consider a sequence of the form (p3) which contains 
j states and (N — j) states \^ 2 ) (he. the number of 
ik which are equal to one is given by j). We denote such 
a sequence by Note that there are b^j sequences 

of this kind. For all those sequences, we find 



Fj = \(*(j)\P A ®P B \*(j))\ 2 



! E 

k-k- <(j+k)<k+ 



Ik 2(N- 



N-j,k\ 



The average fidelity F is given by 



F 



(38) 



(39) 



k~<j<k+ 



where 2~ N is the probability that a certain sequence ap- 
pears, &atj is the number of sequences of the form |W(j)) 
and Fj is given in (p8|). We have that 

1 1 



where — N/2 ± fi/3N@ . In this case, one can also 



bound Fj and finds 



Fj > | £ cfsT-^b 



N-j,k\ , 



k=k- 



where fc ± = c 2 a {N — j) ± /j/3N^ and we have that 
k~ < (j + k) < k + as required. By noting that a bino- 
mial distribution is asymptotically equivalent to a nor- 
mal (Gaussian) distribution, the fidelity Fj Vj can be 
seen to be bounded from below by $(2/xiV /3_1 / 2 ), where 
$(x) = 1/\/2~jt J^ x e v / 2 dy. For our choice of fi,(3, we 
have that Fj — > 1 when N — > oo. In a similar way, one 
shows that also F — > 1 when N — > oo, as after bounding 
Fj as stated above, ( flO| ) also corresponds to a binomial 
distribution centered at j = N/2. This finishes the proof 
of the statements (i-ii). 

In a similar way, one can carry out the analysis for 
a set of L entangled states and an arbitrary probability 
distribution {pt}. In this case, fc ± = N J2t=i Pi^a^^N^ 
and some of the binomial distributions are replaced by 
multinomial distributions. Also in this case, one finds 



that F — > 1 for N — * oo and that the dimension of the 
projector Pa (Pb) is given by the N times the entropy 



of the operator p A ([34 
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